Klein- Tunneling of a Quasirelativistic Bose-Einstein Condensate in an Optical Lattice 
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Optical lattices have proven to be powerful systems for quantum simulations of solid state physics 
effects. Here we report a proof-of-principle experiment simulating effects predicted by relativistic 
wave equations with ultracold atoms in a bichromatic optical lattice that allows for a tailoring of the 
dispersion relation. We observe the analog of Klein-tunneling, the penetration of relativistic particles 
through a potential barrier without the exponential damping that is characteristic for nonrelativistic 
quantum tunneling. Both linear (relativistic) and quadratic (nonrelativistic) dispersion relations are 
investigated, and significant barrier transmission is observed only for the relativistic case. 
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Klein tunneling, a textbook effect in which relativis- 
tic particles penetrate through a potential barrier with- 
out the exponential damping that is characteristic for 
nonrelativistic quantum tunnelling P, 0] , has never been 
observed for elementary particles. In this counterintu- 
itive consequence of relativistic quantum mechanics, a 
strong potential, being repulsive for particles and attrac- 
tive for antiparticles, results in particle- and antiparticle- 
like states aligning in energy across the barrier. There- 
fore, a high transmission probability is expected when 
a potential drop of the order of the particles rest energy 
mc 2 is achieved over the Compton length h/mc. For elec- 
trons, one derives an extremely high required electric field 
strength of s» 10 16 V/cm, which so far has prevented an 
experimental realization on this system. The observation 
of Klein tunnelling has however been reported in solid 
state analogons [3|4g], for example in graphene material. 
In this carbon material owing to a linear, i.e. quasirel- 
ativistic, dispersion relation around the Fermi edge, rel- 
ativistic effects can be very illustratively emulated 
Other systems suitable for the simulation of quasirela- 
tivistic effects are ions with a long-lived two-component 
electronic structure in Paul traps |7|, where the nonrela- 
tivistically forbidden entry into a high potential well has 
been observed Experiments in photonic structures 
and in dark state media 1(J lll have been proposed. 
Ultracold atoms in optical lattices [12] allow for the inves- 
tigation of both linear and nonlinear Hamiltonians due 
to the neutral charge of the atoms, with prospects in- 
cluding the simulation of interacting relativistic quantum 
field theories [13 |. 

Here we report a proof-of-principle quantum simula- 
tion of relativistic wave equation predictions with ultra- 
cold atoms in an optical lattice. Our experiment is based 
on rubidium atoms in a Fourier-synthesized lattice poten- 
tial consisting of an optical standing wave with spatial 
periodicity A/2, where A denotes the laser wavelength, 
and a higher spatial harmonic with A/4 spatial period- 
icity. For a suitable choice of relative phases and am- 
plitudes of the harmonics, the dispersion relation in the 
region between the first and second excited Bloch band 
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FIG. 1. Band structure for an optical lattice potential for a 
relative phase of cp = (left) and <p = n (right). The parame- 
ters for the lattice depths were Vi = 5E r and Vi = 1.6-EV. For 
ip — tv the splitting between the first and the second excited 
band vanishes, and a Dirac-point (marked by the solid blue 
box) is observed. The bands relevant for the experiment are 
shown by solid lines, and on the rescaled energy scale shown 
on the right hand side, zero energy is chosen at the position 
of the band crossing. 



becomes linear, as known for ultrarelativistic particles. 
We experimentally demonstrate both the transmission 
of atoms through a potential barrier for the case of a 
linear dispersion relation, i.e. Klein-tunneling, and the 
usual reflection of atoms by the barrier for the case of a 
quadratic, i.e. Schrodinger-like, dispersion in an excited 
Bloch band. 

The periodic potential used to taylor the dispersion 
of ultracold rubidium atoms is of the form V (z) = 
V 1 /2cos{2kz) + V 2 1 '2 cos (4k z + <p), where k = 2tt/\ is 
the photon wavevector, V\ denotes the potential depth 
of the usual standing wave potential and V 2 that of the 
higher spatial harmonic, generated by the dispersion of 
multiphoton Raman transitions 
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Fig. 1 shows 

the band structure for such a Fourier-synthesized lat- 
tice for V\ = 5E r and V 2 = l.6E r used in our exper- 
iment, where E r = h 2 k 2 /2m denotes the recoil energy, 
for two different values of the relative phase ip between 
lattice harmonics. For the shown parameters, while the 



splitting between the first and the second excited Bloch 
band exhibits a nonzero value for tp = 0, it vanishes for 
a phase of tp = 7r. The critical dependence on the rel- 
ative phase between lattice harmonics is understood in 
terms of the splitting arising from both contributions of 
second order Bragg scattering of the usual lattice and 
of first order Bragg scattering of the higher spatial har- 
monic, where the contributions interfere constructively 
or destructively depending on tp [15j. Of special inter- 
est is the case of destructive interference of the Bragg- 
scattering amplitudes, realized with tp = ir (see Fig. 1 
right), for which at a suitable choice of lattice ampli- 
tudes the dispersion relation near the resulting crossing 
point becomes linear, i.e. relativistic, with an effective 
light speed c e g = 2hk/m ~ l.lcm/s for the used 783nm 
laser wavelength. Both for vanishing and small splittings 
between the bands, we expect to be able to simulate rela- 
tivistic physics, where a variation of the effective atomic 
Compton wavelength X c . e g = h/m e ffc e g — 2c e gh/AE, 
with m e g — AE/2c^ s as the effective mass and AE as 
the size of the splitting, is possible by appropriate choice 
of amplitude and phase values of the lattice harmonics. 
In the limit of AE —> the effective Compton wavelength 
diverges. If we choose the zero point of the energy scale 
to be at the crossing, atomic population in the second ex- 
cited band (above the crossing) corresponds to a particle- 
like excitation, in the below lying band of negative energy 
according to the Stuckelberg-Feynman interpretation to 
a particle-like excitation propagating backwards in time, 
which is equivalent to a temporally forward propagat- 
ing antiparticle excitation Q. Formally, the dynamics 
of atoms in the bichromatic lattice near the crossing be- 
tween the first two excited bands can be described using 
a one-dimensional Dirac-like Hamiltonian (see 16| and 
Supplementary Material): 
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H = m eS cl s a z + c eS qa x + V siow (z) . 
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where a x and a z are Pauli matrices, q — —ih5 z is the 
momentum operator and V s \ ow (z) is an external poten- 
tial varying much slower than the lattice periodicity. 
The two-component Hamiltonian acts on spinors ip = 
(-02 > ^l) , with ipi and ip2 corresponding to course-grain 
atomic wave- functions in the upper and lower bands, 
respectively. Equation [1] becomes exact in the limit 
\q\ <C frk and TO e ffCg ff -C E r . One of the hallmark-effects 
of a relativistic dispersion is Klein-tunneling, which is a 
single-particle effect, so it can be equally well observed 
for bosons (as used in our experiment) and fermions [l7j |. 

In our experiment Klein-tunneling is investigated by 
monitoring the transmission through an external poten- 
tial barrier that stands against the outcoupling of atoms 
from a far detuned optical dipole potential of depth Vq 
by the earth's gravitational acceleration g. The spa- 
tial distribution of the combined potential K;i ow (z) = 

— Vb ex P (—2 (z/wq) 2 ) — mgz is shown in Fig. 2c, where 



(b) 



(c) 



> 




J2. .£ 





position z 

FIG. 2. Klein-tunneling of atoms through a potential barrier, 
(a) Relevant part of the dispersion relation in the lattice for a 
relative phase between lattice harmonics of tp = 0, for which 
the large splitting suppresses a tunneling between bands, (b) 
For ip — 7r the dispersion relation is linear and the two Bloch 
bands touch each other. Shown in the three diagrams is the 
variation of the atomic energy during passage of a potential 
barrier of height Vb- (c) Spatial distribution of the combined 
potential formed by gravitational and dipole trapping poten- 
tial. Atoms can pass the potential well in the case of tp = n 
due to their possibility to drop to below the band crossing 
in the Bloch spectrum when loosing potential energy, see the 
middle graph in (b). 



the height V& of the potential barrier relatively to the 
minimum of the trapping potential can be adjusted by 
variation of the dipole trap beam power. The width of 
the potential barrier is of the same order as the used 
beam diameter 2l}q ~ 46/im. For a typical atomic en- 
ergy of one recoil energy below the maximum of the po- 
tential barrier, the estimated probability for usual non- 
relativistic quantum tunneling is of the order P m — 
exp {—2\Z2mE r z /ft) ~ 10~ 170 , i.e. completely negligi- 
ble. 

The situation however changes when quasirelativistic 
Klein-tunneling occurs, because the tunneling rate for 
this process does not decay exponentially with the width 
of the potential barrier. A quasirelativistic dispersion re- 
lation for ultracold rubidium atoms is induced using the 
Fourier-synthesized optical lattice, and Figs. 2a and 2b 
(left) indicate the relevant part of the atomic dispersion 
relation near the crossing between the first and the sec- 
ond excited Bloch-band for a relative phase of tp = and 
tp = 7r. The atoms are loaded at a quasimomentum q 
well above the crossing region, but when proceeding to- 
wards the potential well on its rising edge loose kinetic 
energy, i.e. their momentum reduces. For a phase tp = 0, 
the splitting between the first and second excited Bloch- 
band is comparatively large (see Fig. 2a). This results 
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in a small effective Compton wavelength, A c e ff ~ 6/zm, 
which is below the length of the rising edge of the bar- 
rier, and we expect no Landau-Zener tunneling into the 
lower band. The dispersion relation for atoms in the up- 
per band then is Schrodinger-like. When the height of 
the potential barrier is larger than q 2 /2m e g, the particle 
cannot pass through the barrier. 

On the other hand, for a relative phase tp = n between 
lattice harmonics, cf. Fig. 2b, the dispersion relation in 
the vicinity of the crossing point between the first and 
the second excited Bloch-band becomes ultrarelativis- 
tic. The effective Compton wavelength becomes larger 
than the widths of the edges of the barrier, and particles 
that approach the potential barrier and loose kinetic en- 
ergy on the rising edge can be accelerated to below the 
crossing point between the second and the first excited 
Bloch-band (the Dirac-point), i.e. to states of negative 
energy of the energy scale shown in Fig. 2c. Correspond- 
ingly, they can surpass higher potential barriers than in 
the nonrelativistic case of Fig. 2a. This corresponds to 
the case of Klein-tunneling. Our experiment simulates 
the conversion of a particle into a spatially backwards 
propagating antiparticlc during the transmission of the 
barrier, whereafter the crossing point of bands again is 
passed on the tailing edge, so that particle-like states are 
again observed beyond the well. As in the case of the 
Klein-tunneling of electrons, this is equivalent to a dou- 
ble Landau-Zener passage between states of positive and 
negative energy respectively [H, Q ■ 

The experiment proceeds by initially producing a Bose- 
Einstein condensate of rubidium atoms in an mj? = 
spin projection of the F = 1 ground state within a dipole 
trapping potential formed by a focused C02-laser beam. 
To prepare atoms in the second excited Bloch-band, we 
initially leave the atoms in ballistic free fall by extin- 
guishing the C02-laser dipole potential until the atoms 
have reached a momentum of 0.9hk by the earths gravi- 
tational force, and then apply a Doppler-sensitive Raman 
pulse transferring atoms to the mp = — 1 spin projection 
and imparting two photon recoil momentum, which in- 
creases the momentum to its desired value of 2.9hk. This 
initial momentum allows to load atoms into the second 
excited Bloch-band of the lattice potential (at the po- 
sition q = 0.9hk), and the C02-laser dipole potential 
required for the shaping of the desired slowly spatially 
varying potential barrier is again activated. The ballistic 
free fall during momentum preparation is about 2/im, i.e. 
well below the focal diameter of the trapping beam. 

Fig. 3a shows typical experimental data for the spa- 
tial atomic distribution, as recorded by absorption imag- 
ing at a time 5ms after preparation, for a relative phase 
tp = and ip = 7r between lattice harmonics respec- 
tively. The height of the potential barrier was V& = 
5E r , while the average initial energy of the particle was 
E = £kj n + w e ff c eff = 3.9_E r , with m e gc^ s of 0.5E r and 
for tp — and tp = it, respectively (see also Fig. la). 
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FIG. 3. (a) Cuts (solid red) through the measured spa- 
tial atomic distribution (insets) for a 5ms experiment time. 
Atoms proceeding from the center of the dipole trapping po- 
tential towards the potential barrier that detains against an 
outcoupling by the earths gravitational field for a relative 
phase between lattice harmonics of (top) tp = and (bottom) 
ip = 7r. The dotted black line indicates the potential Kd OW (z). 
(b) Relative atomic population beyond the potential barrier 
versus phase p. The corresponding effective Compton wave- 
length A Cj0ff = 2c ctt h/AE is shown on the top scale. The 
shown horizontal error bars refer to the latter scale and are 
dominated by the uncertainty in AE, while the estimated un- 
certainty for the phase tp (lower scale) is below the drawing 
size of the dots. The solid line is the result of a numerical in- 
tegration of the relativistic wave equation (see eq. {H). The 
only free fit parameters were amplitude and offset. 



The insets are false-colour shadow images of the atomic 
cloud, with the spot on the left-hand side corresponding 
to atoms near the trap location (i.e. before the barrier) 
and atoms on the right-hand side to particles that have 
transmitted the barrier. The solid red lines in the main 
images are cuts through the center, with the external po- 
tential indicated by dotted black lines. The data shows 
that for a relative phase of tp = 0, almost all atoms re- 
main in front of the barrier, as expected. On the other 
hand, for ip = it, our experimental data show that most 
of the atomic population can be found beyond the poten- 
tial barrier, equivalent to Klein-tunneling in this optical 
lattice system. 

Note that for Klein-tunneling perfect transmission 
through the barrier is expected, when the energy split- 
ting AE would be zero. In our experiment, some 30% of 
the population remains in front of the barrier, which we 
mainly attribute to atoms that did not take part in the 
Doppler-sensitive Raman transfer, and remain trapped in 
the C02-laser beam focus. We have measured the atomic 
population found beyond the barrier for variable values 
of the phase tp between lattice harmonics, to investigate 
the case of intermediate values of the splitting AE be- 
tween Bloch-bands. Corresponding data is displayed by 
the dots in Fig. 3b. The effective Compton wavelength 
becomes sufficiently large for Klein-tunneling only in a 
narrow region near tp = it , corresponding to the ultra- 
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FIG. 4. Fraction of unreflected atoms versus the height of 
the potential barrier Vb for a phase shift (solid red) ip = 
and (dashed blue) if = ir. The solid lines show numerical 
simulations. The grey shaded region corresponds to barrier 
height values where the atomic velocity reaches the edge of 
the Brillouin zone at \q\ = hk. The desired dispersion rela- 
tion is reached only in the white region. The experimental 
parameters are the same as in Fig. 3a except for the barrier 
height. 



relativistic case, while the atoms remain in front of the 
barrier for smaller phase values. The solid line is the 
result of a numerical simulation of the relativistic wave 
equation (see eq.l), in good agreement with the experi- 
mental data. 

A striking prediction for Klein-tunneling is that the 
transmission through the barrier is expected to be inde- 
pendent of the barrier height, an issue in clear contrast 
to the expectations for nonrelativistic quantum mechan- 
ics. For a corresponding measurement in our system a 
momentum resolved time-of-flight measurement was em- 
ployed to allow for a Stern-Gerlach separation of atoms 
that did not take part in the Raman transfer (see Supple- 
mentary Material) . Fig. 4 shows the relative signal of un- 
reflected atoms, corresponding to Klein-tunneled atoms, 
versus the barrier height The grey shaded region 
corresponds to barrier height values where the atomic 
velocity reaches the edge of the lattice Brillouin zone at 
| q | = hk, for which the bottom of the first exited band 
can be reached, i.e. only in the left white region the de- 
sired dispersion relation is achieved, with larger accuracy 
when remaining far from the shaded region. The blue 
squares are data recorded for <p = 7T, corresponding to a 
quasirelativistic dispersion, for which this signal remains 
at a high value within nearly the complete white region, 
illustrating the prediction of Klein-tunneling being inde- 
pendent of the barrier height with good accuracy. On 
the other hand, a pronounced loss of this signal at large 
barrier heights is observed for ip = (red dots), corre- 
sponding to a nonrelativistic dispersion. Quantum tun- 
nelling remains negligible due to the large spatial width 



of the barrier. The finite width of the kinetic atomic 
energy distribution here softens the otherwise expected 
sharp decay for < Vb- The shown solid and dashed 
lines are the result of a numerical simulation, which for 
the relativistic case are in good agreement with the data 
and also qualitatively reproduce the nonrelativistic case. 

To conclude, we report an experiment demonstrating 
the analog of Klein-tunneling, as a proof-of-principle ex- 
periment testing relativistic wave equation predictions, 
with ultracold atoms in a bichromatic optical lattice. By 
tuning the relative phase between lattice harmonics the 
atomic dispersion can be tuned continuously from the 
nonrelativistic to the ultrarelativistic case, though the 
atoms move at a velocity 10 orders of magnitude below 
the speed of light. 

For the future, we expect that ultracold atoms in opti- 
cal lattices allow quantum simulations of a wide range of 
effects of both linear and nonlinear Dirac-dynamics. Per- 
spectives include the verification of theoretical high en- 
ergy physics predictions [3], as chiral confinement, and 
other results of massive Thirring models 13, 3, 20j |. 
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